Local magnetic anisotropy in BaFe 2 As 2 : a polarized inelastic neutron scattering study 
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The anisotropy of the magnetic excitations in BaFe2 As2 was studied by polarized inelastic neutron 
scattering which allows one to separate the components of the magnetic response. Despite the in- 
plane orientation of the static ordered moment we find the in-plane polarized magnons to exhibit a 
larger gap than the out-of-plane polarized ones indicating very strong single-ion anisotropy within 
the layers. It costs more energy to rotate a spin within the orthorhombic a-b plane than rotating it 
perpendicular to the FeAs layers. 
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The recent discovery of the iron-based oxypnictide su- 
perconductorsi has triggered enormous activity explor- 
ing the superconducting pairing mechanism in these ma- 
terials. Density functional theory (DFT) calculations 2,3 
indicate that the conventional electron-phonon coupling 
is by far too weak in the iron pnictides to account for 
the observed high critical temperatures of up to T c — 
55K (Ref.U). Instead, the phase diagrams of many FeAs- 
based families suggest a magnetic pairing mechanism, 
as superconductivity appears close to an antiferromag- 
netic spin-density wave (SDW) phased— Superconduc- 
tivity can even coexist with the SDW order parameter 
which then slightly decreases upon entering the super- 
conducting stated— Furthermore, the onset of super- 
conductivity in doped BaFe2As2 compounds is accompa- 
nied by the appearance of additional magnetic scattering 
at the SDW wave vector*^— 

BaFe2As2 becomes superconducting on hole- 
doping, 7 on electron-doping^ or by the application 
of pressurei 15 i 16 Pure BaFe2As2 exhibits a structural 




FIG. 1: (Color online) (a) Visualization of the crystal and 
magnetic structure of BaFe2As2 in the orthorhombic setting, 
(b) Definition of the reference frame for the neutron polariza- 
tion analysis on the (1 1) reflection as an example (only the 
Fe ions are depicted in the unit cell). 



phase transition from the high-temperature tetragonal 
(14/mmm) to the low-temperature orthorhombic phase 
(Fmmm), see Fig. QJa), at T s -140 K (Ref. [13) 
accompanied by a magnetic phase transition into a 
collincar antiferromagnetic structure with an ordered 
Fe moment of 0.87 [ib (Ref. [l7t ). This orthorhombic 
antiferromagnetic phase exhibits remarkable properties. 
The magnetic structure shown in Fig. 1(b) breaks 
tetragonal symmetry, one may therefore expect some 
orthorhombic splitting in the crystal structure as it 
indeed occurs, but this splitting remains rather small. 
However, the signature of the orthorhombic structure 
in the electronic properties is much stronger than the 
small orthorhombic splitting in the lattice constants 
would suggest. Resistivity measurements on de twinned 
single crystals reveal pronounced anisotropy associ- 
ated already with the precursors of the orthorhombic 
phased The observation of a larger resistivity along the 
shorter axis with ferromagnetic spin coupling is further 
surprising. Also scanning tunnelling microscopy^ and 
optical conductivity analyses 2 - report essential in-plane 
anisotropy. Magnon excitations in the SDW phase have 
been studied by several inelastic neutron scattering 
experiments (unpolarized) and, again, the slope of the 
low energy dispersion is anisotropic The nearest 
neighbor interaction along the a and b directions needed 
to model the magnon dispersion even show opposite 
signs. This remarkable anisotropy in magnetic exchange, 
however, can qualitatively be explained by DFTi2£ The 
pronounced anisotropics in the electronic properties 
inspired a model of orbital polarization 26 which gets 
support from recent ARPES experiments ; ' 28 as well as 
speculations about the role of nematic properties. 
In this work we analyze the magnetic single-ion 
anisotropy which can be deduced from the gap in the 
magnetic excitations finding again remarkably strong 
in-plane anisotropy. It costs more energy to rotate the 
spin within the FeAs layers than aligning it perpendicu- 
lar to them, in contrast to the simple expectation of an 
easy-plane system. 

Large BaFe2As2 single crystals were obtained by a 
self-flux technique. 2 - The pre-reacted precursor materi- 



2 



als FeAs, Fe2As, C02AS and BaAs were mixed leading 
to a Ba(Fe, Co)3,]As3.i composition. This composition 
was used to achieve a homogeneous melt at T = 1463 K. 
The melt was cooled slowly under a temperature gradi- 
ent in a double-wall crucible assembly to obtain large 
and flux- free single crystals of BaFc2As2. Using this 
technique, single crystals with lateral dimensions up to 
25 x 10 mm 2 and thickness up to 1 mm were obtained. 
Two plate-like single crystals with masses of 804 mg and 
200 mg have been coaligned in order to decrease count- 
ing time. The experiment has been carried out at the 
thermal-beam three-axis spectrometer IN20 (ILL) using 
polarizing Heusler (111) crystals as monochromator and 
as analyzer. Longitudinal polarization analysis was per- 
formed with a set of Hclmholtz coils to guide and orient 
the neutron polarization. The longitudinal polarization 
analysis permits separating the nuclear cross-section and 
the individual components of the magnetic cross-section 
into the spin-flip (SF) and non-spin-flip (NSF) channels 
by selecting the initial and the final polarization direction 
of the neutrons. We use the conventional reference frame 
for polarization analysis^ see Fig.QJb), with the x axis 
parallel to the scattering vector Q, y perpendicular to 
x within the scattering plane, and z perpendicular to x 
and y (perpendicular to the scattering plane). The sam- 
ple has been mounted in the [100]/[001] scattering geom- 
etry in orthorhombic notation, but note that the crystal 
is twinned so that part of the crystal is in [010]/[001] 
geometry and signals from both parts superpose. 

Longitudinal polarization analysis adds additional se- 
lection rules to the general law in neutron scattering that 
only magnetic components perpendicular to the scatter- 
ing vector contribute. Only the magnetic components 
perpendicular to the axis of polarization analysis con- 
tribute to the SF scattering whereas the parallel compo- 
nents contribute to the NSF scattering.— Nuclear scat- 
tering is always a NSF process. We measured the NSF 
cross-section for polarization parallel to x, serving as a 
reference for eventual contaminations, and all three SF 
cross-sections keeping constant. Applying a correc- 
tion for the finite flipping ratio of R = 7.4 the magnetic 
cross-sections a y and a z (Eqs. [1] and ^ as well as the 
background in the SF channel (Eq. [3]) can be calculated 
according to: 
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FIG. 2: (Color online) Temperature dependence of the elastic 
magnetic scattering cross section at Q=(l 1) revealing that 
almost the entire magnetic signal stems from a y . The values 
have been obtained by Eqs. [1] and [2] The fitting of a power 
law to the a y data yields a Neel temperature of 137 K. 
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FIG. 3: (Color online) Inelastic constant-Q scans at the mag- 
netic Brillouin-zone centers (1 Z) for different I values and at 
(1 0). (Blue) dots represent the magnetic scattering cross 
section a y , while (red) squares depict a z . The cross sections 
have been obtained by measuring the SF Z and SF y channels, 
respectively, and applying Eqs. [1] and [2] The solid lines are 
spline-interpolated spin-wave theory calculations folded with 
the experimental resolution. 



BGRsf(R) = j^[I(SF v )+I(SF z )} - ^±ll(SF x ). 

(3) 

In Fig. [5] the maximum intensity of the magnetic 
(10 1) reflection is shown as a function of temperature 
for the two channels. A power law has been fitted to the 
Oy data revealing a Tjv of approximately 137 K. The en- 
tire elastic signal is found in a y (a z is essentially zero) in 



perfect agreement with the well-accepted magnetic struc- 
ture shown in Fig. Q] with moments parallel to the a axis. 

Inelastic constant Q-scans were performed at the mag- 
netic zone centers Q=(l I) with i=l,3 and 5, and 
at Q=(l 0) which is a magnetic zone boundary due 
to the finite interlayer coupling (Fig. [3]). The scan at 
Q=(l 0) does not indicate any magnetic scattering 
in agreement with the published magnon dispersion for 
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CaFe 2 As 2 (Refs. MM) and for SrFe 2 As 2 (Refs. HHl) 
indicating energies of the order of 100 meV for this q 
value. Due to the high zone-boundary energies the mag- 
netic signals of the differently oriented twin domains do 
not superimpose in the scans at Q=(l I) with 1=1,3 
and 5, which are zone boundaries for the second domain 
orientation. The magnetic-zone-center scans sense the 
gap of the magnon dispersion arising from anisotropy. 
For a system with uniaxial anisotropy one expects a sin- 
gle gap for the two still degenerate transversally polarized 
magnons. This degeneracy, however, lifts for orthorhom- 
bic local symmetry causing an easy, a medium and a 
hard axis^i We will see that such orthorhombic single- 
ion anisotropy is needed to describe magnetic excitations 
in BaFe 2 As 2 . The cross sections a y and a z unambigu- 
ously show the splitting of the zone-center magnons in 
BaFe 2 As 2 . The z direction of our reference system is 
always parallel to the b direction and therefore a z mea- 
sures the transverse magnon polarized along the b direc- 
tion. The y direction lies always in the a*,c* plane but 
rotates with varying Qi. Since the ordered moment is 
parallel to a, a y contains only the transverse magnon po- 
larized along the c direction (see below for discussion of 
the contribution of the longitudinal magnon). Due to 
the geometry this signal is reduced by a factor 0.92, 0.61, 
and 0.42 for Q;=l,3, and 5, respectively. The qualitative 
analysis of the scans at the magnon zone centers immedi- 
ately shows that the in-plane magnon gap is larger than 
the out-of-plane one. This is remarkable in view of the 
in-plane orientation of the ordered moment^ Starting 
from the tetragonal symmetry at high temperature one 
might describe BaFe 2 As 2 as an easy-plane system similar 
to many layered magnets with K 2 NiF4 (214) structure, 
e.g. K 2 FeF4 (Ref . 1321) . In such a system the two trans- 
verse magnons still split as a higher-order cubic in-plane 
anisotropy adds to the dominating easy-plane anisotropy. 
There is no mass-less Goldstone mode in such system. 
But in this case the energy of the in-plane polarized trans- 
verse magnon is below that of the out-of-plane one^ 
The fact that BaFe 2 As 2 exhibits the opposite behavior, 
the out-of-plane mode lies around 10 meV below the in- 
plane mode near 20 meV, excludes the description of the 
magnetic phase in an adapted tetragonal model, which, 
however, works very well for many 214-materials. The lo- 
cal single-ion anisotropy of the SDW phase in BaFe 2 As 2 
is truly orthorhombic. 3 = This observation is in line with 
the strong electronic signatures^— in this phase which 
contrast with the only weak orthorhombic splitting of the 
lattice. The strong single-ion anisotropy within the FeAs- 
layers gives strong support for pronounced orbital polar- 
ization. Previous unpolarized neutron scattering exper- 
iments on SrFe 2 As 2 , CaFe 2 As 2 and BaFe 2 As 2 (Ref. 1331) 
most likely reported only on the lower gap as the different 
polarizations could not be resolved. 

In order to quantitatively analyze the scattering distri- 
bution we have folded the experimental resolution with 
the scattering function of the magnon contribution us- 
ing the Reslib code.— The dispersion and the structure 
factors of the magnons of the AFe 2 As 2 compounds have 



been calculated in linear spin-wave theory by Yao and 
Carlson 3 - 5 - for the Hamiltonian including nearest-neighbor 
exchange coupling along a and b directions, J a and 
next-nearest neighbor exchange, J 2 , inter- layer coupling, 
J c , and anisotropy terms. A gap at the magnetic zone 
center appears only for a finite single-ion anisotropy term 
J2i ~A(<S?) 2 where, in our case, A is equal 3 ^ to A^ or 
A c and Sf is the component of the spin operator of 
spin i in the direction of the ordered moment. Yao and 
Carlson obtain for the gap at the magnetic zone cen- 
ter: A = 2SVA • (2J a + 4J 2 + A + 2J C ). For the ex- 
change parameters we have used the values recently de- 
termined for SrFe 2 As 2 (Ref. [23h but normalized them to 
the value S—\ -0.87 corresponding to the ordered moment 
in BaFe 2 As 2 , note that the full shape of the dispersion 
is not essential for our analysis as we only sense the low- 
energy part of it in our measurement. With S=^ ■ 0.87, 
S Ji a =30.8 meV, SJ lb =-5 meV, SJ 2 =21.7 meV, 5J C =2.3 
meV, Af,=0.99(2) meV, and A c =0.38(l) meV, we may si- 
multaneously describe the two channels measured at the 
three zone centers, with the anisotropy terms and the 
scale factors as the only free parameters. This descrip- 
tion takes into account the Q-dependence of the mag- 
netic scattering due to the form factor and the geometri- 
cal factor given above. The calculated values have been 
interpolated with a B-spline in order to obtain contin- 
uous curves which are depicted as solid lines in Fig. [3] 
The perfect description of the experimental data vali- 
dates the assumption that the total magnetic scattering 
arises from the transversally polarized magnons. Using 
the obtained anisotropy values we determine the in-plane 
and the out-of-plane gap to amount to 16.4 meV and 10.1 
meV, respectively. 

The single-ion anisotropy parameters obtained in our 
experiment can be compared to the result of relativistic 
DFT calculations including spin-orbit coupling^ the au- 
thors find that only 0.16 meV/Fe are necessary to turn 
the magnetic moment from the easy axis to the b di- 
rection, but that 0.20 meV/Fe are required to turn the 
moment perpendicular to the plane. These calculations 
cannot properly describe the local single-ion anisotropy 
in BaFe 2 As 2 as the sign of the anisotropy does not agree 
with the experimental values. DFT— thus misses some- 
thing in the description of the local symmetry of the SDW 
phase in BaFe 2 As 2 . 

In an itinerant system with non-saturated ordered mo- 
ments, one may expect a longitudinal magnetic excita- 
tion, i.e. a magnon polarized along the ordered moment. 
This mode contributes to the a y channel, but by compar- 
ing the different Q-values studied we can fully rule out 
such an interpretation for the studied energy ranges. The 
longitudinal contribution should increase with increasing 
I component, but the experimental data show the oppo- 
site and are well described by the geometry factors for the 
transverse mode. We can thus excude a longitudinal ex- 
citation below ~35 meV in agreement with DFT studies 
predicting them at significantly higher energiesi 38 i 39 

The temperature dependent measurement of the two 
magnetic cross sections at an energy transfer of 10 meV 
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FIG. 4: (Color online) Temperature dependence of the inelas- 
tic magnetic scattering at Q=(l 1) and an energy trans- 
fer of 10 meV documenting the transition from isotropic to 
anisotropic fluctuations at Tjv. The solid lines serve as a guide 
to the eye. 

illustrates how the isotropic magnetic fluctuations in 
the paramagnetic phase transform into the anisotropic 
magnons below the Neel temperature (Fig. @|. There is 
no difference visible between the two channels at high 
temperature. Upon approaching the magnetic transition 
the intensity increases due to the softening of the mag- 
netic response. Below the transition, the magnetic signal 



is almost completely suppressed in the a z channel where 
the larger gap develops, but there is only moderate re- 
duction in the a y channel, where the energy transfer of 
10 meV is just slightly below the low-temperature maxi- 
mum. 

In conclusion the polarized inelastic neutron scattering 
studies on the magnetic excitations in BaFe2As2 reveal 
a strong in-plane single-ion anisotropy which contrasts 
with the small structural orthorhombic distortion in this 
material. This pronounced local anisotropy corrobo- 
rates the strong electronic signatures of the orthorhombic 
phase and gives further support for an important role of 
orbital degrees of freedom in the iron pnictides. The local 
anisotropy of Fe in the SDW phase apparently is insuf- 
ficiently described in DFT theory even when including 
spin-orbit coupling. 
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